MASTER  COPY 


FOR  REPRODUCTION  PURPOSES 


AD"“  A246  256  *^*^o*t  documentation  page 


2b.  DECLASSIFICATION /OOWNGKAOIMj  SCHEOU 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

Technical  Report  #12 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Dept,  of  Chemistry-SUNYA 


lb.  RESTRICTIVE  MARKING 


3.  DISTRIBUTION /AVAIIA8IUTY  OF  REPORT 

Approved  for  public  release; 
distribution  unlimited. 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6c  ADDRESS  {Ctyv  State,  sni  ZIP  Cotk) 

1400  Washington  Avenue 
Albany,  NY  12222 


6b.  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 
(If  spfilkabit) 

U.  S.  Army  Research  Office 


7b.  ADDRESS  (City,  State,  and  ZIP  Coda) 

P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


Sa.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

U.  S.  Army  Research  Office 


6c  ADDRESS  fCrty,  Stata.  and  ZIP  Coda) 

P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


8b.  OFFICE  SYMBOL  I  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
Of  appOcabla)  I 


to.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 


PROJECT 

NO 


WORK  UNIT 
ACCESSION  NO 


1 1 .  TITLE  (Induda  Saojrity  Classification) 

Diffusion  through  a  membrane  from  a  sessile  droplet-limiting  case  of  the  large  approximation 
for  short  times 


13b.  TIME  COVERED  Il4.  DATE  OF  REPORT  (Taar,  Month,  Day)  |lS.  PAGE  COUNT 

FROM  5/91  TO  8/91  I  1991/8/5 


12.  PERSONAL  AUTHOR(S) 

H.L.  Frisch 


13a.  TYPE  OF  REPORT 

Technical 


16.  SUPPLEMENTARY  NOTATION  „ 

The  view,  opinions  and/or  findings  contained  in  this  report  are  those 
of, the  authqr(s)  and  should  not ,be^ construed  as.  an  official  Department  of  the  Army  position. 


17.  COSATI  COO€S  I  IS.  SUBJECT  TERMS  on  r^vtrso  if  necesMry  #nd  identify  by  bfock  number) 

Diffusion  .droplets,  transport  tbrnugh  membranes- 

ylL  V  ,  .  ,1 


FIELD 

GROUP 

19.  ABSTRACT  {Continue  on  revene  if  neceaery  end  identify  by  block  number) 

^ » 

We  review  briefly  the  singular  perturbation  problem  arising  in  the  approximate 
solution  of  diffusion  from  a  large,  neat,  sessile  agent  droplet  through  a  thin  membrane. 
We  obtain  the  zeroth  order  short  time  solution (which  has  a  definite  thickness  scaling) 
but  this  solution  is  only  valid  for  short  times.  ^  W 

DTIC 

»t.ECT- 

F£Bld19tt|  Ij 


U 


20.  DISTRIBUTION /AVAILABIUTY  OF  ABSTRACT  |21.  ABSTRACT  SECURITY  CLASSIFICATION 

Qunclassifiedainumiteo  □  SAME  as  rpt.  OdticusersI  Unclassified 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL  |22b.  TELEPHONE  (Induda  Area  Coda)  22c.  OFFICE  SYMBOL 


DO  FORM  1473,64  MAR 


B3  APR  tditiofl  may  ba  uMd  until  avbauttad. 
All  otbar  aditiont  are  obtolata. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

UNCLASSIFIED 


Technical  Report  #12 


Diffusion  through  a  membrane  from  a  sessile  droplet-limiting  case 
of  the  large  droplet  approximation  for  short  times. 

by 

H.L.  Frisch 

Department  of  Chemistry 
State  University  of  New  York 
Albany,  NY  12222 


August  5,  1991 


Research  supported  by  the  U.S.  Army  Research  Office  Grant  DAAL  0388K0198 
with  funds  supplied  by  DARPA. 


92  2  12  109 


92-03644 

tl  iSfta  lien  *  ■ 


-2- 


1 .  Introduction 

Recently  I  have  spent  some  effort  on  a  grant  supported  by  the  U.S.  Army 
Materials  and  Mechanics  Research  Center  in  studying  numerically  and  analytically  the 
diffusion  of  agent  from  a  sessile  agent  droplet  (a  spherical  cap)  through  the 
(polymer)  membrane  on  which  it  rests.  The  agent  droplet  can  be  neat  or  polymer 
thickened . ^ ^ ^  Dr.  Nathaniel  Schneider  of  Natick  R.  D.  and  E.  Laboratories  who 
has  asked  us  to  review  these  matters  is  primarily  interested  only  in  the  neat  agent 
problem.  This  diffusion  problem  arises  in  assessing  the  definite  amount  of  agent, 
Q^,  coming  through  a  test  membrane  within  a  definite  time  period,  t^.  Its 
relevance  to  monitoring  and  dosimetry  is  thus  clear  and  needs  no  further  comment. 

^^The  numerical  solutions  make  clear  both  quanlitatively  and  quantitatively  (for 
suitably  chosen  parameters)  the  overall  expected  behavior  and  the  results  are  in 
overall  agreement  with  observations.  While  for  the  polymer  thickened  agent  droplet 
problem  there  is  a  more  extensive  region  of  scaling  with  regard  to  initial  droplet 
radius,  R  ,  and  membrane  thickness,  L,  this  is  not  the  case  for  the  neat 

Q  .  , 

i  1-3 

droplet.  Specifically  we  showed-  that  simple  scaling  relations  apply  only 
for  a  rather  limited  parameter  ranges:  a)  The  more  ubiquitous  (under  field 
conditions)  case  where  the  agent  is  deposited  as  an  aerosol  of  tiny  droplets  which 


spread  on  the  membrane,  for  j;^ich  R^/L  «  l.-^the  small  droplet  case^  This  case  is 
essentially  completely  analyzed  in  references  (1)  and  (A)  and  a  quasi  steady  state 

analytical  theory  is  given  in  reference  (2)  for  the  case  that  the  agent  solubility  - 

is  sufficiently  large.  I  am  given  to  understand  by  Dr.  Schneider  that  the  scaling 


$ 


and  applicability  of  this  case  is  sufficiently  well  understood  that  no  further 

comment  is  necessary  here,  (b)  The  more  special  case,  we  call  the  laree  droplet 

case,  arises  When  R  /L  »  1.  Note  that  the  condition  used  for  approximate 

o 


AvaUablilty  Cod^ 
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solutlon  of  this  problem  that  the  droplet  is  smaller  than  L  cannot  be  correct  for 
all  time.^  The  radius  of  the  sessile  droplet  decreases  continuously  with  time  and 
will  ultimately  become  equat  to  or  smaller  than  L  (as  the  agent  diffuses  through  the 
membrane)  when  the  special  scaling  for  the  limiting  ease  must  fail.  In  the  next 
section  we  review  the  discussion  of  this  special  limiting  case.  It  has  been 
originally  treated  by  us  in  reference  (1)  and  (3).  We  shall  treat  it  here  more 
extensively  in  the  next  section  showing  that  it  arises  from  a  singular  perturbation 
of  the  original  boundary  value  problem.  In  the  last  section  we  deal  with  the 
extension  of  this  problem  when  evaporation  from  the  sessile  droplet  surface  is  not 
negligible. 

2.  The  large  droplet  case  -  negligible  surface  evaporation  from  the  droplet. 

We  consider  the  moving  boundary  diffusion  problem  initially  described  in 

references  (1)  -  (4).  We  use  a  slightly  modified  notation.  Ml  physical  entities 

with  dimensions,  other  than  certain  scaling  variables  such  as  R^,  L,  etc.  will  be 

written  with  a  carat  over  the  letter  denoting  the  entity.  We  will  subsequently 

introduce  suitably  scaled,  dimensionless  variables  denoted  by  letters  without  carats. 

The  concentration  c  of  agent  in  the  membrane  at  location  r,  z  at  time  t 

1-4 

satisfies  the  diffusion  boundary  value  problem  : 

Oc/at)  =  D  {(32c/ar2)  +  1/?  (ac/ar)  +  (a2S/az2)}  d) 
inO<z<L  ,  r>0,t>0 

with  initial  and  boundary  conditions 

c(r,  r,  t  =  0)  =  0,  (2) 

c(r,  0,  t)  =  c  ,  r  <  R(t)  (3) 

o  ~ 

with  c  the  solubility  of  agent,  i.e.  the  equilibrium  concentration  of  agent  in  a 
o 

membrane  immersed  in  neat  liquid  agent  at  ambient  fixed  temperature  and  pressure  and 

A  A 

R(t)  the  time  dependent  sessile  droplet  radius.  The  sessile  droplet  is  a  spherical 
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cap  with  fixed  contact  angle  9  and  volume  1/3  f(e)R  (t), 

_2 

f(9)  a  «sin9(2-i-cos6)/(l-fcos9)  .  Furthermore, 

Oc/3z)2_q  =  0  ,  r  >  R(t) 


Ur,  L,  t)  a  0.  (5) 

A  A  A 

The  moving  boundary  of  the  spherical  cap  r  =  R(t)  is  determined  by  the  mass  balance 

A  A 

A  A  A  A 

pf(e)R2  (dR/dt)  a  2irDfdr  r  (3c/dz)$^0  •  R(t=0)  =  Rq  ,  (6) 

0 

with  p  the  density  of  neat  agent.  We  now  introduce  dimensionless  variables: 

A  A 

r  a  r/R^  ,  z  a  z/L, 

t  a  Dt/L^  ,  R(t)  a  R(t)/R  , 

o 

c(r,z,t)  a  c<r,z,t)/c^, 

6  a  2»CQ/pf(e)  ,  (7) 

and  a  dimensionless  parameter  of  smallness, 

c  a  L/R  «  1, 
o 

into  this  boundary  value  problem. 

One  finds 

(3c/3t)  a  (3^c/3z^)  +  [(3^c/3r^)  +  1/r  3c/3r)]  (8) 

in0<z<l,  r>0,  t>0; 

c(r,  z,  0)  a  0  ,  (9) 


c(r,0,t)  a  1  for  r  <  R(t), 


(3c/3z)  =•  0  for  r  >  R(t)  ; 

ZaO 

c(r,l,t)  a  0 
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R(t) 

r2  (dR/dt)  =.  c  B  J  dr  r  (ac/az)^-©  ,  (13) 

0 

R(0)  -  1. 

The  natural  way  to  seek  an  approximate  solution  of  (8)  -  (13)  is  to  consider  the 
singular  perturbation  in  which  the  zeroth  order  solution  for  the  concentration 
c^  3  c^(z,t)  satisfies  (8)  with  neglect  of  the  term  in  the  square  bracket  on  the 
r.h.s.  of  (8),  i.e. 

c  3  c„(z,t)  +  0(c)  (14) 

o 

R(t)  3  R^(t)  +  cR^(t)  +  0(c) 
where 

ac  /at  3  a^c  /az^  ,  o  <  z  <  i  ,  t  >  o; 
o  o 

c^(z,0)  3  0  .  (15) 

c^(0,t)  3  1  , 

c^d.t)  -  0; 
and  using  (14)  in  (13) 

dR  /dt  3  0  ,  R  3  1  ,  (16) 

o  o 

Rqsi 

(dRi/dt)  3  Bj  dr  r  (aco/3z)o  =  B/2  (aCo/dz)o  ,  (17) 

0 

Rl(0)  3  0 

The  zeroth  order  solution  of  the  boundary  value  problem  (cf.  (14)  and  (15))  has  been 
given  by  us  in  reference  (1)  and  also  in  reference  (3)  where  it  is  given  by  the 
special  case  described  by  Eq.  [7]  on  p.  284  which  in  our  notation  reads: 

Co(z.t)  3  t  {erfc  (2n  +  z/lfl.)  -  erfc  ( (2(n+l)-z/2  ft)  1 }  (18) 

n3i 

t 

R(t)  3  1  +  (cB/2)  J  dt’  I(aCo(z,t’)/dz)l  2.0  ,  (19) 

0 
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The  integrals  given  above  are  tabulated  in  standard  tables. 

To  obtain  a  perturbation  solution  Which  holds  for  longer  times  one  must 
introduce  "stretched”  coordinates,  and  find  the  "inner"  and  "outer"  solution  of  this 
singular  perturbation.  We  shall  not  carry  out  this  rather  dull  but  lengthy  standard 
problem  in  applied  mathematics. 

Instead  we  point  out  that  the  zeroth  order  solution  provides  a  limiting  case, 
valid  for  very  short  times,  when  we  further  can  replace  the  infinite  series  in  (18) 
by  its  leading  term  for  which 

(3Co/3z)2=0~  -I(irt)-l-/2  .  (20) 

Substituting  this  into  (17)  yields 

-  /1.x  or/i./  xl/2  ,.,/.  1  xl/2  -1/t 

R^(t)  =  -a[(t/ir)  +  {2(t/ir)  e 

2  erfc  (l//t)}]  (21) 

or 

R(t)  ~  1-  eB((t/ir)^^^  +  {2(t/ir) -  2  erfcd/ /t) } I . 

In  zeroth  approximation  the  flux  through  themembrane  at  z  =  1,  (which  in  this 
approximation  is  exactly  the  same  as  the  area  averaged  flux  over  the  membrane)  can 
be  obtained  from  (18).  We  shall  again  only  retain  the  leading  order  term  in  the 
summation  (18)  to  calculate  the  zeroth  order  flux,  which  we  write  in  dimensional 


form. 


F(t)  -  -D  (3c  /3z)/'  , 

o  z=L  ^ 

o  /«/  i.xl^2  -if/ADt 
».  2c  (D/»t)  c 
o 


The  absence  of  R  in  (22)  clearly  shows  that  (22)  can  only  apply  for  short 
o 


enough  times  even  if  c  <  <  1.  The  use  of  (22)  to  calculate  the  reduced, 

/'  2 

dimensionless  break  through  time,  t.  D/L  ,  requires  that 

b 

A  2 

t.  D/L  <  T, 

D  L 


(23a) 
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where  t  is  the  time  needed  for  the  droplet  radius  to  decrease  by  L,  i.e. 

R(t  )  a  1-e.  In  zeroth  order  approximation  the  time  t  is  defined  by  the 
L  L 

solution  of  the  transcendental  equation 

B  -2  erfc(l/  V^)}] 

=  1.  (23b) 

The  special  scaling  in  (22)  also  noted  by  Angelopoulos,  Schneider  and  Meldon 

subsequently  to  us,  in  an  unpublished  report  applies  only  well  for  sufficiently 

short  times,  even  for  R  /L  >  1. 

o 

3.  The  zeroth  order  solution  of  the  large  droplet  case  for  non-negligible 
evaporation. 

As  shown  in  reference  (4)  if  evaporation  is  not  negligible  (4)  has  to  be 
replaced  by  (dc/dz)A_Q  +  h  [c(r,0,t)-U)  =  0  ,  r  >  R(t)  (24) 

and  (b)  by 


A  A 


0 


-R2f(e)yp 

A 

with  h  the  evaporation  coefficient  characterizing  the  flat  polymer-vapor  interface, 

^  A 

U  the  equilibrium  concentration  of  agent  in  the  vapor  above  the  polymer,  y  the 
(per  unit  surface  area)  evaporation  coefficient  of  agent  from  a  neat  agent  surface 
at  ambient  conditions.  In  the  zeroth  order  approximation  we  neglect  all  diffusion 
for  r  >  R(t)  so  that  the  solution  c^(z,t)  is  unaffected  by  (24)  and  is  still  given 
by  (18).  Denoting  by  y  the  appropriately  scaled,  dimensionless  y  defined  by 
y  =  Ly/D 


(26) 
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we  can  write  the  dimensionless  equation  (25)  as 

R 

r2  (dR/dt)  »  c{fl/dr  r  Oc/az)2=0  <27) 

0 

Prom  (27)  and  (14),  Ro(t)  =  1,  but 
1 

(dR^/dt)  =  Bf  dr  r  Oco/az)2,o  -2yRi  (28) 

0 

with  Ri(0)  s  0.  The  solution  of  (28)  is 


Rl(t)  »  fljdf  e-2Y(t-f)  (3co(r,z,f)/az)z=o  (29) 


which  can  be  evaluated  easily  for  short  times  using  (20),  the  resulting  integrals 
are  all  tabulated  . 
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